Probing a critical length scale at the glass transition 
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We give evidence of a clear structural signature of the glass transition, in terms of a static correla- 
tion length with the same dependence on the system size which is typical of critical phenomena. Our 
approach is to introduce an external, static perturbation to extract the structural information from 
the system's response. In particular, we consider the transformation behavior of the local minima 
of the underlying potential energy landscape (inherent structures), under a static deformation. The 
finite-size scaling analysis of our numerical results indicate that the correlation length diverges at 
a temperature Tc, below the temperatures where the system can be equilibrated. Our numerical 
results are consistent with random first order theory, which predicts such a divergence with a critical 
exponent v — 2/3 at the Kauzmann temperature, where the extrapolated configurational entropy 
vanishes. 

PACS numbers: 64.70.kj,61.43.Fs,64.70.Q-,05.20.Jj 
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Is there a static, structural origin of the dramatic slow- 
ing down of dynamics in supercooled liquids? The phys- 
ical mechanism of the dramatic slowing down and corre- 
sponding increase in viscosity upon approaching the glass 
transition has been vehemently debated for more than 40 
years [H [5] . In critical phenomena, the slowing down of 
the system's dynamics is intimately related to a nearby 
phase transition, where a diverging length scale typically 
reflects the onset of long-range correlation. Several re- 
cent theories of the glass transition are built in this spirit 
[SHS] and are supported from results on dielectric suscep- 
tibility [7] and specific heat [SI IS] ■ The few experimen- 
tal studies, however, show only a moderate cooperativity 
length near the glass transition of some four to five parti- 
cle diameters [10] . Specifically designed simulations [11] 
have recently indicated a growing amorphous order in the 
range of four particle diameters. The results confirm, to 
some extent, the predictions of random first order theory 
(RFOT) [3l, f6] . However, no evidence for critical behav- 
ior with a diverging length scale could be found there. A 
growing correlation length is also indicated by a very re- 
cent numerical study of different local quantities |12j and 
it is interpreted as the hint of a critical behavior, whose 
evidence, nevertheless, is still very elusive. The hetero- 
geneous dynamics of glass forming liquids allow one to 
define a dynamical correlation length [TJ [T3l - fT6] which is 
indeed growing as the glass transition is approached, but 
it remains unclear how, or whether at all, this dynam- 
ical length has a static, structural origin as is the case 
in phase transitions. A promising, alternative approach 
to supercooled liquids is based on inherent structures, 
i.e. the local minima of the underlying potential energy 
landscape [T7]. Although several qualitative changes in 
the system dynamics can be related to changes in their 
inherent structures [THl [IH], no length scale as signifi- 
cantly increasing as the dynamical correlation length has 
emerged so far [20] . 

Following a nonequilibrium thermodynamic theory of 
glasses [21] , here we use small, static deformations to per- 
turb the inherent structure configurations of supercooled 



liquids approaching the glass transition [55]. With this 
procedure, large correlated regions emerge at low temper- 
atures and allow us to detect a static correlation length 
that shows critical behavior upon approaching the glass 
transition. 

Methods and numerical simulation. - We employ a 
binary Lennard- Jones mixture, which is an established 
model for fragile glass formers [52] (see Appendix IA| . 
An ensemble of well equilibrated configurations are pre- 
pared at constant density by slowly cooling statistically 
independent samples from high temperatures down to 
the supercooled regime. We study different system sizes 
ranging from N = 2000 up to iV == 64000 particles. The 
inherent structure X^ = {r'^} corresponding to the ac- 
tual configuration X — {tj} is obtained by locally mini- 
mizing the system's potential energy by a conjugate gra- 
dient method. The potential energy landscape and its 
local minima (the inherent structures) have been stud- 
ied intensively in recent years, unraveling a number of 
remarkable relations between inherent structure proper- 
ties and the system's behavior [51H55I . Unlike these pre- 
vious studies, we investigate the relation between two 
inherent structure configurations, X'^ and X'^'^, where 
the latter is the inherent structure corresponding to the 
affinely deformed configuration X'^ = {r^}, rf = E • r^, 
where the deformation is represented by the matrix E. 
From these configurations, we define the mismatch vec- 
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E • r^, which give the nonaffine displace- 
ments between two inherent structure configurations. In 
particular, we consider static shear deformations with 
small amplitude 7, E = 1-1- 76162, where 6^ are Carte- 
sian unit vectors. In Ref. [55], we found, as the tem- 
perature is lowered towards the glass transition, char- 
acteristic changes which are strongly reminiscent of the 
systems long time dynamics: small average mismatch 
lengths which are correlated over large distances. 

Results. - In order to study correlations in the direc- 
tion of the mismatch field, we define the coarse-grained 
nonaffine displacement field D;,(r) = ^ Dj{b)S{r — r^) 
of the inherent structure configuration X^, where Dj(6) 







FIG. 1: (color online) . Coarse-grained nonaffine displacement 
field D6(r) defined in Eq. Q for high (T = 1.0, left panel) 
and low (T = 0.4, right panel) temperature. The coarse- 
graining length was chosen as 6 = 2.0. The color code illus- 
trates surfaces of constant |Dp-values. For better visibility, 
only particles with |Djp > 0.5 are shown. 



are the coarse-grained mismatch orientations obtained by 
averaging d^- 



dj7|dj| over a sphere of radius b 



i^^(^) = ^7'EdU''0 



,dq\ 
jk) 
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Here, r ■? is the distance between particles j and k in the 
inherent structure of deformed configuration X'^'^, Nj — 
Sfc X6('',fe) the number of neighbors of particles j within 
a distance b, and Xb{i^) = 1 if ?' < 6 and zero elsewhere. 
For sufficiently small deformations, the mismatch field 
is approximately independent of strain amplitude 7 in 
a statistical sense and the use of normalized mismatch 
vectors d' suppresses the strong temperature dependence 
of the average mismatch lengths [22] . 

In Fig. [1] Df,(r) of a typical configuration is shown 
at high and low temperature. The large domains where 
|Djp > 0.5 at low temperature clearly indicate large 
regions where the response to the quasi-static external 
deformation is strongly correlated, in marked contrast to 
the high temperature behavior. 

These observations can be made more quantitative 
by studying the distribution of coarse-grained mismatch 
lengths /ib(D^). For small b, only a single particle con- 
tributes to ^ and /ib(D^) -^ (5(D^ — 1). For increasing 6, 
Dj includes more and more particles and therefore, with 
the decay of correlations, /if, accumulates more weight at 
small values of D^ until /i(,(D^) — > <5(D^) when b becomes 
large. The transition between mostly ordered (peak of 
ht, near 1) and disordered (peak near 0) regions happens 
at a value of b which strongly increases with decreasing 
temperature (see Fig. I2|. 

The coarse-graining length b is our means to probe the 
extent of correlations in the mismatch field on lengths 
up to 6. Therefore, from the histogram /ib(D^) we de- 
fine a static correlation length ^b as the coarse-graining 
distance b for which the peak-location has moved from 
one to 1/e. Our results are robust and the conclusions 
presented here hold also for other definitions of S^b, see 
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FIG. 2: Histogram of mismatch lengths ^i,(D^) at high (a: 
T — 1.0) and low (b: T = 0.4) temperature for different 
values of the coarse-graining length b. At high temperatures, 
the typical value of D^ jumps from an ordered (D'^ ~ 1) to 
a disordered one (D^ « 0) when 6 becomes slightly larger 
than the particle diameter. At low temperatures, instead, 
this transition happens at much larger values of 6. 



Appendix |B] Figure [3] shows the correlation length ^b 
as a function of temperature. In the high temperature 
liquid phase we observe a low value of ^s of the order of 
the particles diameter a, weakly dependent on temper- 
ature. In the supercooled regime, however, ^b increases 
considerably. For the lowest temperatures at which we 
can equilibrate the system, T « 0.40, we find ^s « 4ct 
for the largest systems studied here. 

We also observe in Fig. [3] a systematic finite-size effect 
oi^B{T) in the supercooled regime, where the correlation 
length grows with the linear size L of the system. We 
investigate this dependence in more detail because, to- 
gether with the significant increase of (,b, this behavior is 
typical of critical phenomena |30| , where the correlation 
length ^ associated with the order parameter character- 
izing the transition diverges at the critical temperature 
Tc as £, ^ {T — Tc)^'^ ■ Therefore, we perform the same 
finite-size analysis which in critical phenomena is used to 
extract the critical behavior from data obtained in a finite 
system. If the increase of ^b upon lowering the tempera- 
ture is indeed due to a diverging static correlation length 
f underlying the glass transition, then ^b should diverge 
as well, and in the infinite system ^^(r) ^ t^'', where 
t — {T — Tc) is the distance from the critical point Tc 
and p the critical exponent. On the basis of the scaling 
hypothesis for critical phenomena [30! , the corresponding 
quantity S,b,l{T) in a finite system of linear size L should 
follow the behavior 
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where Q^^ {x) is a universal scaling function and v is 
the critical exponent associated to ^. Therefore, we 
plot ^b.l{T)L~p/'^ as a function of the scaling variable 
(T - TjLi/'^ in Fig. [3} We identify the critical temper- 
ature Tc with the Kauzmann temperature T/^, where a 
phase transition due to an entropy crisis should be lo- 
cated according to RFOT and related theories [3l |6] . For 
the system used here, Tfj- has been determined numeri- 
cally as Tk ~ 0.30 [HI [3T]. Fixing this value for T^ we 
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FIG. 3: (color online). Top: Static correlation length ^s of 
coarse-grained mismatch vectors as function of temperature 
for different system sizes. Bottom: Data collapse of £,b,l{T) 
according to Eq. (pi. 



observe that the numerical data for different tempera- 
tures T and system sizes L collapse onto a single master 
curve in agreement with Eq. ([2]). We find the best data 
collapse for the critical exponents p « 0.9 ±0.1, and 

V K, 0.65 ± 0.1, which is the case shown in Fig. |3] The 
value of V is very close to u = 2/3, which is indeed pre- 
dicted by RFOT for three-dimensional systems [Sj and to 

V w 0.69 deduced from light scattering experiments on o- 
terphenyl [10] . The value v = 2/3 is also compatible with 
recent simulations |12| . where, however, no proper finite- 
size scaling was performed. From the quality of data 
collapse and the intersection of ^b{T)L^''^^ for different 
system sizes [30], we estimate that the critical tempera- 
ture lies in the range 0.25 < T^ < 0.4. Using the finite- 
size analysis in this range of temperatures we end up with 
the estimates for the critical exponents p w 0.8 ±0.2 and 

V K, 0.7 ± 0.15. The precision of these results is lim- 
ited by the increase of ^b being too weak for the smallest 
system sizes, which is probably also the reason why previ- 
ous studies could extract only small, roughly temperature 
independent length scales from inherent structures [24] . 
Moreover, the finite-size analysis in the vicinity of the 
critical point is made extremely hard by the long equili- 
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FIG. 4: (a) The structural relaxation time t is plotted ver- 
sus static correlation length ^s for the system with A'' = 
64000 particles. The solid line is the best fit to the data 
and shows an exponential dependence of the form r ~ 
exp [^'(^s/fcflT)^'^'']. (b) The static correlation length is 
shown as a function of the configurational entropy 5conf ■ The 
solid line is the best fit to the data and shows the power-law 
relation fs ~ (TSco„f)~'°^ 



bration times. In spite of these limitations, the critical 
region is apparently large enough to be felt in the acces- 
sible temperature regime. 

Within RFOT, the increasing relaxation time r 
is linked to the free energy barrier of nucleating a 
new structure of linear size ^m in the liquid, t ^ 
exp[B(T5conf)-'/(''-'')] - e^^[A{U/kBTf] ^mM- 
The configurational entropy per particle is denoted by 
Sconf , d is the spatial dimension and 9 = d—l/v the ex- 
ponent of the interface free energy cost of a nucleus. This 
gives us the possibility to further check the consistency 
of our data with the RFOT scenario. In Fig. |4|-a, we 
plot the relaxation time r versus ^s/fc^r on a double- 
logarithmic scale. The relaxation times are large-A^ val- 
ues of T taken from Ref. [5D], which have been deter- 
mined therein from the final relaxation of the incoherent 
scattering function. The fitted line shows that the form 
T - exp [A{i/kBTY] - exp [A^^B/kBTf"'"] indeed de- 
scribes our numerical data quite well. The value of the 
exponent 9 we obtain is 6* sa 1.6 ± 0.3. If we identify our 
correlation length ^ with the mosaic length scale ^m in 
RFOT, this value of 9 we find is in agreement with the 
prediction 9 = d/2 of RFOT [3 . It is worth noting that 
our analysis is free of inconsistencies reported in Ref. [50] 
when looking at length scales from dynamical heterogene- 
ity. We can also directly test the connection between 
configurational entropy Scant and correlation length ^, 
which according to RFOT reads ^ - (TS'conf)^^/^''"^^ 
In Fig. |4]-b, ^b is plotted as a function of TSconf on a 
double-logarithmic scale. The data for the latter are 
taken from Ref. [32]. Our results indicate that the re- 

P — p 

lation ^B ^ i" ^ {TSconi)"^''^''^ is indeed vahd for the 
largest system we studied. For the exponent, we find 
{d — 9)~^ « 0.74 ± 0.14 which again is in agreement with 
the value 2/(i predicted by RFOT and which in addition 
implies the validity of the Vogel-Fulcher-Tamman expres- 



sion r - exp [DT/{T - Tk)] ^■ 

Discussion. - Bringing together concepts from inher- 
ent structure formahsm of supercooled hquids and from 
deformations of amorphous sohds offers a new perspec- 
tive on the long-standing glass problem. In this spirit, we 
have employed a new method in order to elucidate grow- 
ing static correlations in supercooled liquids. We found 
direct evidence for the critical behavior of the static cor- 
relation length close to the glass transition, which has 
been speculated by several recent theories [SHS]. In par- 
ticular, our estimates of the transition temperature, crit- 
ical exponents, as well as the link with relaxational dy- 
namics are in quantitative agreement with the predic- 
tions from RFOT [3] . However further studies are needed 
in order to determine the numerical values more pre- 
cisely. Our results have important consequences for flow- 
induced rearrangements in colloidal glasses [33] and elas- 
tic response of amorphous solids [29] , where similar co- 
operative behavior has been observed recently. 

Our study has been performed at constant density, 
where the critical temperature Tc extracted from ^b is 
apparently below the lowest temperature at which we 
have been able to equilibrate the system. However, the 
critical point would not necessarily be located at the den- 
sity chosen here and only a systematic investigation in the 
(T, P)-plane will shed further light on its existence, na- 
ture, and location. The long-standing problem of a true 
thermodynamic phase transition underlying the slowing 
down of the dynamics at the glass transition can now be 
addressed from a different perspective, knowing that a 
direct structural signature exists. 



Appendix A: Details of Numerical Simulations 

In this study, we employ the well-known Kob- Andersen 
binary Lennard- Jones mixture 23J, which is a classical 
model for a glass forming liquid. Both types of particles 
{A and B) have the same unit mass and all particles 
interact with each other via a Lennard- Jones potential, 



Vabir) ^ Aeat,[{cJat,/rY'' ~ [aab/rf 



(Al) 



a, 6 e {A, B}, where r denotes the distance between par- 
ticles. For computational convenience, the interactions 
are cut off at r^ = 2.5(Tab and shifted so that Vab{rc) = 0. 
The interaction energies and diameters for the different 
species are given by aAA — c, ctab = O.Sa, asB = 0.88cr, 
<^AA = e, ^AB — l-5e, and sbb — 0.5e. The compo- 
sition of A and B particles is chosen as 80:20 and the 
particle density p = N/V is p = 1.2 /a^. Simulations 
are performed in the A^T^T-ensemble, for a cubic simula- 
tion box subject to periodic boundary conditions. Stan- 
dard molecular dynamics simulations are performed with 
the simulation package LAMMPS [M], using the veloc- 
ity Verlet integrator with a time step increasing with 
temperature from At = 0.001 tref to 0.005 irefj where 
trof = (mcTQ/e)^'^ is the reference Lennard- Jones time. 
An ensemble of statistically independent configurations 



are generated at high temperatures and cooled slowly, 
by coupling the system to a Nose-Hoover thermostat 
with prescribed cooling protocol. The lowest tempera- 
ture configurations have been equilibrated for 5 x 10^ 
MD steps. We use dimensionless temperature T mea- 
sured in energy units e. For the finite-size analysis we 
use system consisting of TV = 2000, 4000, 8000, 32000 
and 64000 particles, which correspond to system sizes of 
L = ^1/3 ^ 11.8,14.9,18.8,29.84, and 37.6 ct, respec- 
tively. 

We obtain the inherent structure configurations X"^ 
from the equilibrated MD configurations X by lo- 
cally minimizing the total potential energy Vtot = 
J2i<j^i^ij)^ using the conjugate gradient method. We 
use the Polak-Ribiere version of the algorithm as imple- 
mented in LAMMPS [34]. The minimization is stopped 
when the potential energy change is less than a tolerance 
value 10~^e. We verified that the results are insensitive 
to a further decrease of the tolerance level. 



Appendix B: Definition of correlation length 

The correlation length S^b can be defined from the dis- 
tribution of coarse-grained mismatch vectors in slightly 
different ways. In the manuscript, we have identi- 
fied £,B with that value of 6, where the peak location 
of hb has moved from 1 to 1/e. Alternatively, we 
have also considered that from the mean-squared lengths 
Jyhb{y)dy — B{bY one can obtain the so-called coarse- 
graining function B{b), defined as the average mismatch 
length when smeared out over a sphere of radius 6, 
B{b) = (Dj(&)2)i/2^ as in Ref. [gg-. By construction, 
5(0) = 1 and _B — > when the coarse-graining length 
6 spans the whole system, since there is no overall mis- 
match, X^i '^7" ~ 0- I^or high temperatures, B{b) de- 
cays so rapidly that we find it difficult to extract reli- 
able values for the correlation length. For relatively low 
temperatures, we observe an exponential decay B{b) ex 
exp [—6/(2^^)], with a well-defined correlation length ^^. 
This regime is similar to the case of amorphous solids [29] . 
We have introduced a factor 2 in the definition of ^^ for 
a closer similarity to the choice in Ref. [5^ , where coarse- 
graining over cubic boxes of side lengths b was performed. 
For a sharply peaked distribution hf){y), ^b — Cb- I^ ^he 
general case, they will give slightly different values but 
we have verified that the temperature dependence of ^g 
in the low temperature regime is fully consistent with the 
results extracted from Fig. 3. Finally, we have also con- 
sidered the length scale ^^ — 61/3, where 61/3 is the value 
of b where one-third of the particles have D^ > 1/2 and 
two-thirds D^ < 1/2. We expect ^g to be close to ^b but 
somewhat smoother and less noisy since integration can 
remove ambiguities in locating the peak position. Also in 
this case the temperature dependence is fully consistent 
with the one of ^b ■ 

Figure SI shows the correlation lengths ^b,S,'b'^b ^*^" 
cording to the above definitions as a function of tempera- 
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FIG. 5: (color online). Correlation lengths CsiCsifs defined 
in the text as a function of temperature T. The system con- 
tains N = 64000 particles. 



ture. The system size considered is TV = 64000 particles. 
All three definitions give consistent results. 
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